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ABSTRACT 

We calculate the type la supernova rate for different star formation histories 
in galaxies by adopting the most popular and recent progenitor models. We 
show that the timescale for the maximum in the type la supernova rate, which 
corresponds also to time of the maximum enrichment, is not unique but is a 
strong function of the adopted stellar lifetimes, initial mass function and star 
formation rate. This timescale varies from ~ 40 — 50 Myr for an instantaneous 
starburst to ~ 0.3 Gyr for a typical elliptical galaxy to ~ 4.0 — 5.0 Gyr for a disk 
of a spiral Galaxy like the Milky Way. We also show that the typical timescale 
of 1 Gyr, often quoted as the typical timescale for the type la supernovae, is 
just the time at which, in the solar neighbourhood, the Fe production from 
supernovae la starts to become important and not the time at which SNe la 
start to explode. As a cosequence of this, a change in slope in the [O/Fe] ratio 
is expected in correspondance of this timescale. We conclude that the suggested 
lack of supernovae la at low metallicities produces results at variance with the 
observed [O/Fe] vs. [Fe/H] relation in the solar region. We also compute the 
supernova la rates for different galaxies as a function of redshift and predict an 
extended maximum between redshift z ~ 3.6 and z ~ 1.6 for elliptical galaxies, 
and two maxima, one at z ~ 3 and the other aX z ~ 1, for spiral galaxies, 
under the assumption that galaxies start forming stars at Zf ~ 5 and Om = 0.3, 
tt A = 0.7. 

Subject headings: supernovae: general- Stars: binaries-nucleosynthesis 

1. INTRODUCTION 

The supernovae of type la are crucial in understanding a number of astrophysical prob- 
lems of primary importance, such as the supernova (SN) progenitors, the determinations 
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of cosmological constants, the chemical enrichment of galaxies and the thermal history of 
the interstellar (ISM) and intracluster (ICM) medium. In the past year there has been a 
great deal of work concerning the different roles of different types of SNe on the evolution of 
galaxies. 

Tinsley (1980) first discussed the possibility of iron being mostly produced in SNe of 
type I (the distinction in type la and lb was not common at that time) and in the following 
years type la supernovae started to be considered as important iron contributors in the 
chemical evolution of galaxies (Greggio and Renzini 1983a,b; Matteucci and Greggio 1986). 
In particular, the observed [O/Fe] versus [Fe/H] relation was interpreted as due to the 
different roles played by SNe of type II (explosion of single massive stars) and type la in the 
production of oxygen and iron. Supernovae of type II are, in fact, responsible for the bulk of 
oxygen and for part of iron production whereas type la are responsible for a large fraction 
of iron production, the exact value depending on the assumed initial mass function. 

At the same time models for the progenitors of SNe of type la were developed (Iben and 
Tutukov 1984; 1985). Such models suggest that these SNe may originate from white dwarfs 
in binary systems, either systems with two white dwarfs or systems with one white dwarf 
and a red (super) giant, as originally proposed by Whelan and Iben (1973). In parallel, 
models of explosive nucleosynthesis for such systems (mostly the C-deflagration/detonation 
of a C-0 white dwarf of Chandrasekhar mass) were also presented: the most famous of these 
models being the W7 model of Nomoto et al. (1984). Calculations of the type la SN rates 
were performed on the basis of these models (Greggio and Renzini 1983b; Matteucci and 
Greggio 1986; Matteucci and Tornambe 1987; Tornambe and Matteucci 1987; Matteucci and 
Tornambe 1988) both for our galaxy and for elliptical galaxies. In particular Tornambe and 
Matteucci (1987) and Matteucci and Tornambe (1988) showed that the predicted type la 
rates for our Galaxy and elliptical galaxies are very different because of the differences in 
their histories of star formation. They showed that the maximum of the type la SN rate in 
ellipticals occurs much earlier than in the solar neighbourhood. The same result was shown 
more recently by Matteucci (1994) who adopted a more sophisticated model, including dark 
matter, for the evolution of elliptical galaxies. This result has strong implications on the 
definition of the typical timescale for the chemical enrichment by type la SNe. Yoshii et 
al. (1996) pointed out that the typical timescale for the enrichment from SN la in the solar 
vicinity can be identified with the time at which the [O/Fe] ratio starts changing its slope. 
Independently of the assumed model for progenitors of type la SNe, these authors identified 
that this timescale for the solar neighbourhood is tj a ~ 1.5 Gyr, in agreement with previous 
estimates by Matteucci and Frangois (1992). By coincidence, the value of [Fe/H] at which 
the [O/Fe] ratio changes slope coincides with the transition between halo and disk ([Fe/H]~ 
-1.0). Therefore, this timescale can also be taken as an indicator for the duration of the halo- 
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thick disk phase. However, tj a cannot be universal, as it is often claimed in the literature, 
but it depends strongly on the assumed star formation history. As a consequence of that, 
we expect that the [a/Fe] vs. [Fe/H] relationships to be different in galaxies with different 
histories of star formation (Matteucci 1991). 

We will compute here the typical timescale for the maximum enrichment from type la 
SNe, tj a , in galaxies with different histories of star formation (starburst, elliptical and spiral 
galaxies), defined as the time of the maximum in the type la SN rate. In particular, in the 
framework of different star formation histories, we will explore different progenitor models 
for type la SNe, such as the modification proposed by Greggio (1996) to the classical rate 
computed by Greggio and Renzini (1983b) and the more recent model proposed by Hachisu 
et al. (1996; 1999). 

In section 2 we recall the most popular models for the progenitors of type la SNe and 
define ti a . In section 3 we compute the rates of SNe of type la in galaxies of different 
morphological type. In section 4 we discuss the effects of different type la SN rates on 
galactic chemical evolution. Finally, in section 5 some conclusions are drawn. 

2. TYPE IA SNe, THEIR PROGENITORS, THEIR TIMESCALES 

2.1. Supernovae la 

The main observational constraints about SNela are the lack of hydrogen lines, the 
presence of Sill lines, together with some other intermediate-mass elements (Call, Sn, Ol), 
with a wide range of velocities (~ 10 4 — 3 x 10 4 Km s _1 ), which dominate the photospheric 
spectra near the time of the maximum light. In the later nebular phase, the spectra are 
dominated by iron features, thus indicating that the inner layers consist of iron-peak elements 
(~ 0.1 - 1M ). 

Most of the observed SNIa lie in a narrow range of parameters (spectra, light curves, 
absolute magnitude peak) and only ~ 10 % of the observational sample are conspicuously 
peculiar (mainly subluminous). 

In the light of these constraints, it is widely accepted that SNela originates from the 
thermonuclear disruption of a white dwarf (WD) which accretes material from a companion 
star in close binary systems, but there are still a number of uncertainties about the nature of 
the companion (main sequence star, red giant or another WD), the nature of the white dwarf 
(C-0 WD, He WD or O-Ne-Mg WD), the mass reached by the accretor at the explosion 
(Chandrasekhar or sub-Chandrasekhar) and the explosion mechanism itself (deflagration, 
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delayed-detonation or detonation). However, there is a general consensus in assuming that 
a C-deflagration in a C-0 WD of Chandrasekhar mass best represents the characteristics of 
the majority of the SNe la. 

2.2. Progenitors 

We recall here the most common models for the progenitors of type la SNe proposed 
insofar: 

• The merging of two C-0 WDs, due to gravitational wave radiation, which reach the 
Chandrasekhar mass and explode by C-deflagration (Iben and Tutukov 1984; 1985). 
This is known as double-degenerate (DD) scenario. 

• The C-deflagration of a Chandrasekhar mass (~ 1.4M ) C-0 WD after accretion 
from a non-degenerate companion (Whelan and Iben 1973; Munari and Renzini 1992; 
Kenyon et al. 1993). This model is known as the single-degenerate (SD) one. The 
main problem with this scenario is the narrow range of permitted values of the mass 
accretion rate in order to obtain a stable accretion, instead of an unstable accretion 
with a consequent nova explosion and mass loss. In this case, in fact, the WD never 
achieves the Chandrasekhar mass. In particular, Nomoto, Thielemann & Yokoi (1984) 
found that a central carbon-deflagration of a WD results for a high accretion rate 
(M > 4 • 1CT 8 M yr" 1 ) from the secondary to the primary star (the WD). They found 
that ~ 0.6 — 0.7 M Q of Fe plus traces of elements from C to Si are produced in the 
deflagration, well reproducing the observed spectra. 

• A sub- Chandrasekhar C-0 WD exploding by He-detonation induced by accretion of 
He-rich material from a He star companion (Limongi and Tornambe 1991). 

• A recent model by Hachisu et al. (1996; 1999) is based on the classical scenario of 
Whelan and Iben (1973) (namely C-deflagration in a WD reaching the Chandrasekhar 
mass after accreting material from a star which fills its Roche lobe), but they find 
an important metallicity effect. When the accretion process begins, the primary star 
(WD) develops an optically thick wind which helps in stabilizing the mass transfer 
process. When the metallicity is low ([Fe/H] < — 1), the stellar wind is too weak and 
the explosion cannot occur. This model is appealing since it overcomes the difficulty 
of achieving a stable accretion as in the classic SD scenario discussed above. However, 
this model still needs to be tested in several astrophysical contexts, as we will see in 
the next sections. 
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It is worth noting that in the last few years the DD scenario has lost some credibility 
with respect to the SD scenario, mostly because of the negative results of observational 
searches for very close binary systems made of massive enough white dwarfs (Bragaglia et 
al. 1990). Therefore, in this paper we will concentrate on the SD model. 



2.3. Timescales 

In the framework of the SD scenario, the explosion timescales for different progenitor 
models are given by the lifetime of the secondary star. In particular: 

• In the formulation of the type la rate by Greggio and Renzini (1983b) (hereafter 
GR83), based on the Whelan and Iben (1973) model (hereafter WI73) the explosion 
times correspond to the lifetimes of stars in the mass range 0.8 — 8M Q . In fact, the 
maximum initial mass which leads to the formation of a C-0 WD is ~ 8M Q , although 
stellar models with overshooting predict a lower value (e.g. Marigo et al. 1996), which 
means that the first system, made of two 8M stars, explodes after ~ 3 • 10 7 years from 
the beginning of star formation. The minimum mass of the binary system is assumed 
to be 3 M , to ensure that the WD and the companion are large enough to allow the 
C-0 WD with the minimum possible mass (~ 0.5M Q , since stars with masses below 
this limit can only become He-WDs) to reach the Chandrasekhar mass after accretion. 
The smallest possible secondary mass is O.8M and therefore the maximum explosion 
time is the age of the universe. This ensures that this model is able to predict a present 
time SN la rate for those galaxies where star formation must have stopped several Gyr 
ago, such as ellipticals. 

• Greggio (1996) (hereafter G96) revised the computation of the SN la rate in the frame- 
work of the SD model and suggested a more detailed criterium for the formation of a 
system which can eventually explode as type la SN. In particular, the explosion will 
occur when: 

M WD +eM 2 , e >M Ch , (1) 

where M 2 e is the envelope mass of the evolving secondary and e is the accretion 
efficiency, namely the mass fraction of the envelope which accretes over the WD. The 
maximum value for the efficiency is e — 1. The mass Mwd is the mass of the white 
dwarf: we have assumed three different Mj n — Mwd relationships, from Renzini and 
Voli (1981) (RV81), from Iben (1991) (191) and from Marigo et al. (1996; 1998) (MBC). 
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• Kobayashi et al. (1998, 2000) (hereafter K98 and KTN) adopted Hachisu et al. (1996; 
1999) model and considered two possible progenitor systems: either a WD plus a red 
giant (RG) star or a WD plus a main sequence star (MS). In both cases the mass of 
the primary star is defined in the range 3.0 — 8M , whereas the secondary masses are 



(1.8M© < M MS ^ 2.6M ) and (O.9M < M RG < 1.5M ) for the MS and the RG stars, 



respectively. In this case, the explosion times are given again simply by the lifetimes 
of the secondary stars. If we compare this model with the GR83 model, we see that 
the explosion timescales here are much longer, since the most massive secondary is a 
2.6M star with a lifetime of ~ 3.3 • 10 8 years, instead of a 8M with a lifetime of 
~ 3.0 ■ 10 7 years. In addition, the systems able to give rise to SNe la may form only 
after the gas has reached a metallicity of [Fe/H] = —1.0 dex (K98) or [Fe/H] = —1.1 
dex (KTN) , and this introduces an additional delay which has important consequences 
on galactic chemical evolution as we will see in the next sections. 



We adopt as definition of the typical timescale for the type la SN maximum enrichment 
the time tj a at which the SNIa rate reaches a maximum. As we will see in the next sections 
this maximum depends upon the assumed type la SN progenitor model, namely the lifetimes 
of the progenitors and the initial mass function (IMF), and the star formation rate (SFR). 
The combination of all these parameters creates different shapes and different maxima in 
the type la SN rate versus time relationship. 



In this section we describe the computation of the SN la rate according to different 
histories of star formation: a) an instantaneous burst (one stellar generation), b) constant 
star formation rate, c) a more realistic star formation rate derived in the framework of 
chemical evolution models. 

• a) The calculation of the SNIa rate following a burst of star formation of negligible 
duration, in the framework of the SD scenario and the WI73 model, can be expressed 



2.4. Definition of tj, 



a 



3. 



THEORETICAL SN IA RATES 



as: 
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where M B = Mi + M 2 is the total mass of the binary system, M B)inf and M BjSup are 
the minimum and maximum masses for the binary systems contributing at the time 
t. The maximum value that M B can assume is called Mbm and the minimum M Bm . 
These values (maximum and minimum mass of the binary system able to produce 
a SNIa explosion) are model-dependent. In particular, GR83 considered that only 
stars with M < 8M Q could develop a degenerate C-0 core, thus obtaining an upper 
limit M BM = 16M© for the mass of the binary system. The adopted lower limit is 
(M Bm = 3M ), as discussed previously. The extremes of the integral (2), at a fixed 
time t, are: 

M B>in f = max{2M 2 {t), M Bm ) (3) 

M B , sup = ^M BM + M 2 (t), (4) 
where M BjSup = M BM when M 2 {t) = 8M . 

We define \x = M 2 /M B as the mass fraction of the secondary and f(/j,) is the distribution 
function of this ratio. Statistical studies (e.g. Tutukov & Yungelson 1980) indicates 
that mass ratios close to one are preferred, so the formula: 

/(/i) = 2 1+ ^(l + 7 )^, (5) 

is commonly adopted (GR83), with 7=2 as a parameter. However, it is worth men- 
tioning that in the literature there are different suggestions about the value of 7. For 
example, Duquennoy and Mayor (1991) found a value of 7=-0.35, the value adopted 
in K98. Therefore, we will explore in this paper more values for the parameter 7 (see 
Table 1). 

In order to obtain the SNIa rate, we need also the function M 2 {t). To this purpose we 
adopt the inverse of the formula: 



t(M) = io[ 1 - 338 -v /L79 ~ - 2232 ( 7 - 764 - lo s( M ))]/ - 1116 yr , (6) 
(Padovani & Matteucci 1993). 

K98 adopted stellar lifetimes depending on the initial stellar metallicities from Kodama 
(1997): these timescales are systematically longer than ours especially in the domain 
of low and intermediate mass stars. 

The function <f>(M B ) is the IMF: 
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0(M B ) oc M^ L+X) (7) 

where x — 1.35 is the Salpeter (1955) index and the IMF is defined in the mass interval 
0.1-100 M Q . In the paper we have used both the Salpeter and the Scalo(1986) IMF. 
Finally, the constant A, already defined in GR83, is the fraction of binary systems, 
in the IMF, with the right characteristics to become type la SNe and is fixed by 
reproducing the present time observed type la SN rate. 

For the G96 model the expression of the SN la rate is the same as eq. (2) whereas the 
lower limit for the binary system M B;in f at the time t is: 

M B ,i„f = max(2 M 2 (t), M lmin + M 2 (t)), (8) 

where M lmin is the minimum mass the primary should have to give rise to a WD 
whose mass satisfies eq. (1). We consider also a Sub-Chandrasekhar model in which 
the explosion occurred when Mwd > O.6M and e ■ M 2fi > 0.15 M . For this model 
we assume e = 0.5. 

The calculation of the SNIa rate with KTN prescriptions is a little more complicated, 
because we have different restrictions for the primary and the secondary star: the 
secondary lies in a narrow range of masses (O.9M < M 2 < 1.5M for the WD+RG 
system and 1.8M < M 2 < 2.6M for the WD+MS system). Then, K98 and KTN 
consider only WDs with masses 0.7 M < M WD < 1.2 M , which corresponds to an 
interval for the primary star 3 M < Mi < 8M . In this case the maximum and 
minimum values allowed for the progenitor binary systems should be computed by 
taking into account the above conditions, namely, Mbm = 9.5M and M Bm = 3.9M 
for the case WD + RG and M BM = 1O.6M and M Bm = 4.8M for the case WD + 
MS. For the mass M B)in f we have: 

M B)inf = max(2 M 2 , M Bm , M lmm + M 2 ) (9) 
and for the mass M BjSup : 

M BjSup = min(M BM , M lmax + M 2 ) (10) 
where M lmax = 8M and M lmin = 3M . 

In Figure 1 we show the SN la rate computed in the case of an instantaneous burst for 
the GR83, G96 and KTN (without the metallicity effect) prescriptions. As one can see, 
in the first two cases the rate first increases, reaches a maximum at ti a ~ 5 • 10 7 years 
and then declines. This behaviour is due to the competition between the increasing 
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number of stars with decreasing mass (IMF) and the increasing stellar lifetime with 
decreasing mass. The maximum in the G96 cases is slightly anticipated relative to 
GR83, owing to the fact that a maximum mass producing a C-0 WD of 8.8 M & was 
adopted. The other differences between the G96 (e = 1, 7 = 0.5) and GR83 (7 = 0.5) 
cases are due to condition (1) and the different Mj — Mwd relations adopted. On the 
other hand, the G96 (e < 1) cases predict quite different SN la rates. In particular, 
as the efficiency of mass transfer decreases, the number of systems able to produce a 
type la SN is drastically reduced because of condition (1). As a consequence of this, 
the case with e = 0.2 predicts that type la SN rate following an instantaneous burst 
goes to zero already after a time smaller than 1 Gyr. 

The rate predicted by the KTN model shows instead a completely different behaviour, 
with a discontinuity due to jump from 1.8 M to 1.5 M Q in the mass of the secondary 
passing from the system WD + MS to WD + RG. In figure 1 the rates are normalized 
to their own maximum value as in GR83, whereas in figure 2 we plot the absolute SN 
rates normalized to reproduce a short starburst like those occurring in blue compact 
galaxies. In particular, we assumed A=0.006, a value which can reproduce the features 
of IZwl8 (see Recchi et al. 2001). 

• b) The single-burst approximation could be valid only for starburst galaxies whereas 
for elliptical and spiral galaxies we must consider a star formation extended in time. 
In the case in which the SFR, ^(t), is constant in time, we obtain: 

R la (t)=A^> 0(M B ) / f(fi)dfidM B , (11) 

J A4"B,inf J /l min 

Where \i is defined in the range /i m i n — 0.5, with: 

M 2 M B -0.5M BM 
= max( — , — ) (12) 

for GR83 and G96, whereas for KTN we define: 

M 2 (M B -M lmax ) M 2 
u min = maxi— — , — , — — — J (13) 

^ M B M B M B ' v ' 

and: 

. ,M2max {Mb — Mimin) 

/W = ™n(— , Mb } (14) 

The equation (11) is useful to calculate the SNIa rate in the GR83 and G96 cases 
either for a burst with constant star formation rate or for the solar neighbourhood 
where the SFR has not varied much in time, whereas in all the other cases a constant 
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SFR represent a poor approximation. For a correct use of the assumption of K98 and 
KTN we must consider the metal enrichment of the ISM (we cannot produce SNIa 
progenitors until [Fe/H] = — 1.0 (K98) or —1.1 (KTN). So, in order to calculate the 
SNIa rate, either with constant or variable SFR, we need a chemical evolution model, 
as shown in c). 

• c) The case with a more realistic star formation rate first increasing then reaching a 
maxiumum and decreasing, as predicted by successful models for the chemical evolution 
of the solar neighbourhood (e.g. Chiappini et al. 1997), can be written as: 



as shown first by Matteucci & Greggio (1986). The star formation rate in this case 
has to be evaluated at the time (t — tm 2 ), with tm 2 being the clock for the explosion. 
In figure 3 we show the type la SN rates normalized to their maxima as obtained by 
eq. (15) when the SFR of the Chiappini et al. (1997) model is adopted. Such a SFR 
is proportional to a power k = 1.5 of the surface gas density and to a power h = 0.5 of 
the total surface mass density. Such a formulation of the SFR takes into account the 
feedback mechanism between stars and gas regulating star formation and is supported 
by observations (e.g. Dopita and Ryder 1994). We adopted two chemical evolution 
models: i) the one of Matteucci and Frangois (1992) which is similar to that adopted 
by K98 and KTN. In this model the halo and disk form out of a unique infall episode 
with a time scale for the formation of the disk at the solar neighbourhood of r = 4 
Gyr. ii) the model of Chiappini et al. (1997) which is known as two infall model where 
the halo and thick disk are assumed to form on a short timescale (~ 1 Gyr) out of 
a first episode of infall of primordial gas, whereas the thin-disk is assumed to have 
formed on a longer timescale (~ 8 Gyr at the solar circle). In this model, the existence 
of a threshold in the gas density to regulate the star formation (Kennicutt 1989; 1998) 
is also adopted. As a consequence of this, the star formation rate goes to zero every 
time that the gas density decreases below the threshold (~ 7M Q pc~ 2 ). This model 
reproduces the majority of the features of the solar vicinity and the whole disk and we 
consider it as the best model. In all the models we adopted the SFR of Chiappini et 
al. (1997) as described before, but only in one case we assume the threshold in the gas 
density. 

The model parameters are described in Table 1, where the second column indicates if 
the chemical evolution model for the solar neighbourhood has one or two infall episodes. 
The third column shows the efficiency of star formation v (i.e. the constant in the SFR, 
expressed in Gyr -1 ) and the fourth the adopted timescale (in Gyr) for the formation 
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of the disk in the solar region. The fifth column contains the adopted values for the 
parameter 7 in the distribution function of the mass ratios in binary systems. The sixth 
column indicates if a threshold in the gas density for the SFR has been adopted or not. 
The seventh column shows the assumed prescriptions for the type la SN progenitors 
in the various models. The eighth, ninth and tenth columns contain the values for A, 
IMF and M{ — Mwd relation, respectively. Finally, in the eleventh column is shown 
the present time SN la rate as predicted by each model. The yields adopted in all 
the models are the same: RV81 for low and intermediate mass stars and Woosley and 
Weaver (1995) for massive stars. 

As one can see, in figure 3 the SNIa rate as predicted with the KTN prescriptions 
including the metallicity effect (model Ml) shows a maximum very late (t Ia ~ 8 Gyr), 
as it should be expected given the nature of the assumed progenitors. It is worth noting 
that the model parameters in Ml are the same adopted in the chemical evolution 
model described in KTN. The model with one infall but the prescriptions of GR83 
(M3) predicts a maximum at around tj a = 4 — 5 Gyr, whereas the two-infall model 
(M2) predictions are a bit more complicated showing two maxima, the first one at 
ti a = 1.5 Gyr, due to the first infall episode and the associated star formation, and 
the second one between 4 and 5 Gyr due to the SN la systems born at the beginning 
of the thin-disk formation. The one infall model with G96 prescriptions (M4) presents 
a maximum at t Ia = 3 Gyr. In figure 3 is shown also the SN la rate as predicted 
by Matteucci (1994) by adopting the GR83 prescription, for an elliptical galaxy with 
10 11 Mq of luminous mass. In this case, the SFR is very efficient (~ 10 — 15 times 
more than in the solar neighbourhood) and lasts only for ~ 0.4 Gyr. Therefore, the 
maximum of the SN la rate occurs already at tj a = 0.3 Gyr. In figure 4 the absolute 
type la SN rates, computed by means of the models of Table 1, are presented. In order 
to compute the absolute rates we had to assume a specific value for the constant A, 
which represents the fraction of binary systems in the IMF developing a type la SN 
explosion. In principle, the parameter A should be fixed by reproducing the observed 
SN la rate but it should also assume reasonable values. In our best model (M2), which 
adopts a Scalo (1986) IMF with two slopes (x 1 = 1.35 for M < 2M Q and x 2 = 1.7 
for M > 2M ), the adopted value is A = 0.05 and the predicted present time type 
la SN rate is in good agreement with the observed one (~ 0.18 SN 100 yr~ l , with 
Hq = 60 Km sec" 1 Mpc' 1 , Cappellaro et al. 1999), as shown in Table 1. For Ml, 
with the prescriptions of KTN, we adopted, as in their paper, the Salpeter IMF with 
Ams — 0.05 and A RG = 0.02 (where A MS refers to the systems WD+MS and A RG 
refers to the systems WD+RG), in order to be able to compare our chemical results 
with theirs. The predicted present time SN la rate in this case is too low. We could 
obtain a reasonable value for the SN la rate with the K98 and KTN prescriptions 
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only by assuming A = 0.3, which is a very high and perhaps unrealistic value for such 
a fraction. In fact, a fraction of 30 % of binary systems giving rise to type la SNe 
would imply a total fraction of interacting binaries in the mass range producing the 
SNe larger than 50 % (De Donder and Vanbeveren 2001). The same problem arises 
with the G96 prescriptions, as we can see in figure 4 where the rates of type la SN, as 
predicted by models M4 and M5, are shown. In fact, models M4 and M5 predict too 
low type la SN rates at the present time unless we assume A = 0.3. The differences 
between M4 and M5 are due only to the different M, L — M WD relations adopted in the 
two cases. 



4. THE EFFECTS OF THE TYPE IA SN RATE ON CHEMICAL 

EVOLUTION 

In this section we discuss the impact of the different SN la rates we have discussed up 
to now on the [O/Fe] versus [Fe/H] relationship, which is often used in the literature to infer 
the timescale for SN la enrichment. In figure 5 we show the data from Gratton et al. (2000) 
compared with different model predictions (in all the models the normalization is made with 
the observed solar values of Anders and Grevesse, 1989) including those taken from the paper 
of K98. These latter seem to be the best in reproducing the observed trend but this result 
is not clear. In fact, the predictions of our Ml (see figure 5), which is the copy of KTN 
model, show a very different behaviour which is, in our opinion, more understandable: in 
particular, it does not show any change in slope in the [O/Fe] ratio over the whole range of 
[Fe/H], since the first systems giving rise to type la SNe form only when the gas reaches the 
metallicity [Fe/H] =-1.1, and therefore they will contribute to the chemical enrichment much 
later than the time corresponding to [Fe/H] =-1.1 dex and even [Fe/H]=-1.0 dex, as indicated 
by the maximum of the SN la rate occurring at 8 Gyr. We do not understand how K98, 
who assumed [Fe/H] =-1.0 dex as the limiting metallicity for the formation of the type la SN 
progenitors, can find a change in slope right at this precise metallicity. We conclude that the 
type la SN rate suggested by K98 and KTN, besides giving too low values for the present 
time rate, unless an enormous amount of binary systems is assumed, cannot reproduce the 
[O/Fe] versus [Fe/H] relationship. This relationship, in fact, is characterized observationally 
by a marked change in slope occurring at around [Fe/H] =-1.0 dex, which corresponds to the 
transition between halo and disk metallicities. In our best model the change in slope occurs 
roughly at this point and it corresponds to a galactic age of ~ 1 Gyr. This timescale, in turn, 
corresponds to the first maximum reached by the type la SN rate in the two-infall model and 
can be considered as the typical timescale for the SN la enrichment to become important. It 
is worth noting that in our best model the trend of the [O/Fe] ratio for [Fe/H] < 1.0 dex is 
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not flat but it shows a slight slope since the very first type la SNe start occurring already at 
30 Myr from the beginning of star formation (GR83 model), as discussed also by Chiappini 
et al. (1999). 

In conclusion, the best prescriptions for progenitor models of type la SNe in order to 
reproduce the observed type la SN rate and the chemical evolution of the Galaxy is still the 
GR83 model. 



5. CONCLUSIONS 

In this paper we have explored various SN la progenitor models and the resulting SN 
rates for different histories of star formation and we have defined the typical timescale for 
the maximum chemical enrichment from type la SNe in different galaxies. We have studied 
the effects of the SD scenario, as originally proposed by WI73, with the recipes of GR83 
and G96 as well as the more recent SD model proposed by K98 and KTN. Then we have 
calculated the chemical evolution of the solar neighbourhood in the different cases by means 
of chemical evolution models. Our results can be summarized as follows: 

• The best prescriptions to obtain type la SN rates in agreement with the observations 
seem to be those of GR83. If one adopts the more realistic approach developed by 
G96 to calculate the rate of explosion of systems made of a C-0 WD plus a red giant 
star, then a large fraction (~ 30 %) of interacting binary systems giving rise to type 
la SNe is required to reproduce the observed rate and the solar Fe abundance. This 
corresponds to an overall interacting binary frequency > 50% in the mass range of the 
primary star. The same is true for the model K98 and KTN, which presents also some 
other problems concerning galactic chemical evolution. 

• In particular, the models of K98 and KTN do not produce satisfactory results for the 
[O/Fe] vs. [Fe/H] relation in the solar neighbourhood since it predicts that binary 
systems able to form type la SNe can form very late, after the gas in the Galaxy 
has reached [Fe/H] =-1.1- -1.0 dex. This fact clearly prevents the change in slope 
in the [O/Fe] ratio right at [Fe/H] ~ -1.0, as observed in the data. The predicted 
metallicity effect seems also at variance with recent observations of Damped Lyman- 
a systems (e.g. Pettini et al. 1999) indicating almost solar [cc/Fe] ratios in objects 
with metallicities lower than [Fe/H]=-1.0 dex, thus suggesting that type la SNe have 
already polluted the interstellar medium at low metallicities. However, KTN suggested 
a possible explanation for this effect. In particular, they attributed the almost solar 
ratios at low metallicities as due to inhomogeneous chemical evolution in those early 
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phases. In our opinion this solution looks rather "ad hoc" and is not justifiable for 
metallicities [Fe/H] > —3.0 dex, since no large spread is observed for these metallicities 
in the abundance ratios of halo stars (Argast et al. 2000) 

• The typical timescale for the maximum chemical enrichment by SN la, tj a , can be 
defined as the time at which the SN la rate reaches a maximum. This timescale 
depends upon the adopted SN progenitor model, the stellar lifetimes, the IMF and the 
SFR. Since galaxies of different morphological type are likely to have different histories 
of star formation, the typical type la SN timescale are necessarily different from galaxy 
to galaxy. 

• We have studied different cases in the framework of GR83 model for SN la progenitors: 
in an instantaneous starburst, which can approximate the situation of blue compact 
galaxies, (t Ia ) B c ~ 40 — 50 Myr, in a spiral galaxy like the Milky Way, where the SFR 
has not varied much in time but it has been more ore less continuous, (t/ a ) $ P ~ 4 — 5 
Gyr. A minor maximum appears at ~ 1 Gyr in the two-infall model, where the 
halo is explicitely taken into account and which better fits the properties of the solar 
neighbourhood, corresponding to the bulk of type la SNe originating from the star 
formation in the halo. In an elliptical galaxy, where the SFR was quite high but lasted 
for a relatively short time (~ 0.2 — 0.4 Gyr), (t/ a )s ~ 0.3 Gyr. 

• The knee observed in the [a/Fe] vs. [Fe/H] relations in the solar neighbourhood stars 
corresponds to a timescale t^nee = 1 Gyr. This is the time at which the Fe production 
from type la SNe starts to become important, due to the systems formed during the 
halo phase, and is not a universal value but peculiar to the solar region. 

• The type la SN rates per unit luminosity were plotted as function of redshift (see figure 
6), in a cosmology where we assume that all the galaxies formed at Zf = 5, Qu = 0.3 
and Q\ = 0.7. The Hubble constant is assumed to be H = 60 Km s~ x Mpc~ x . The 
photometric model we adopted is that of Jimenez et al. (1998) As one can see in figure 
6, where the rates for a typical elliptical and a typical spiral galaxy, expressed in SNu, 
are shown, a maximum at a redshift z between 3 and 4 and extending until z ~ 1.6 
should be expected for the elliptical galaxies, whereas we predict two peaks for spiral 
galaxies (two-infall model) the first at z ~ 3 during the halo phase and the second 
at z ~ 1, during the disk phase. In the same figure are indicated the present time 
values for the type la rate in ellipticals and spirals as well as the type la rate derived 
at redshift z ~ 0.4 by Pain et al. (1996). The observed values have been scaled to 
our adopted H . The agreement between the observed and predicted rates seems quite 
good (especially for the present time rates), although it would be more appropriate to 
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compute the cosmic SN rate (i.e. weighted over the different morphological types) in 
order to compare it with the high redshift observations. 
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Table 1: Model Parameters and Predicted SN la rates 



Model 


Infall 


V 


T 


7 


Threshold 


SNela 


A 


IMF 


Mi - M WD 


Rates (lOOyr- 1 ) 


Ml 


1 


0.37 


5 


-0.35 


No 


KTN 


0.05 WDH 


-MS 


Salpeter* 


RV81 


0.025 
















0.02 WDH 


-RG 








M2 


2 


2; 1 


8 


2 


Yes 


C97 


0.05 




Scalo 


RV81 


0.172 


M3 


1 


0.5 


4 


2 


No 


C97 


0.05 




Scalo 


RV81 


0.193 


M4 


1 


0.5 


4 


0.5 


No 


G96 


0.05 




Scalo 


191 


0.039 


M5 


1 


0.5 


4 


0.5 


No 


G96 


0.05 




Scalo 


MBC 


0.023 



* defined in the range 0.05-50 M ( 
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8 9 10 

log t(yr) 

Fig. I. — SNIa rates vs. time for an instantaneous burst of star formation. These rates 
are normalized to their respective maximum values. The solid line is the GR83 model. The 
dotted lines are the G96 models with different values of e (0.2, 0.35, 0.5 and 1 respectively), 
the long-dashed line is the Sub-Ch G96 model (e = 0.5) and the short-dashed lines are the 
K98 model without the metallicity effect. In all the models we have assumed a Salpeter IMF 
and 7 = 2.0, except in the K98 model where we assumed 7 = —0.35 in order to reproduce 
exactly their model. A model with G96 prescriptions and 7 = 0.5 is also shown (dot-dashed 
line). The M; — M WD relationship for the solid line is from RV81 whereas for the G96 model 
is from 191. 
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Fig. 2. — Same as Fig. 1, but with the absolute values of SNIa rates, in units of SNe yr . 
The assumed value of A=0.006 is from Recchi et al. (2001) and is chosen to fit the properties 
of the galaxy IZwl8. 
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Fig. 3. — SNIa rates normalized to their maxima and computed by means of a star formation 
rate depending on the gas density in the framework of the chemical evolution models shown 
in Table 1. The continuous line corresponds to Ml, the short-dashed line to M3 and the line 
with the small oscillations to M2, where a threshold in the gas density for the star formation 
is assumed. The long-dashed line corresponds to model M4. The short-dashed-long-dashed 
line is the rate predicted for an elliptical galaxy by Matteucci (1994). 
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Fig. 4. — The absolute SN rates in units of SNe yr^ 1 computed for the models described in 
Table 1, where all the assumed model parameters are shown. The symbols are the same as 
in fig. 3. We report here also the prediction of Model 5 (dot-dashed line). The G96 models 
(Mod. 4 and Mod. 5) are computed by adopting the e — 1. 
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Fig. 5. — The observed and predicted [O/Fe] vs. [Fe/H] relationships. The predictions 
are from the models of Table 1. The abundances are all normalised to the observed solar 
abundances (Anders and Grevesse, 1989). The dashed-dotted line is the model computed by 
K98, their figure 3, that we show for comparison. The data points (filled squares) are from 
Gratton et al. (2000). 
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Fig. 6. — Predicted type la SN rates expressed in SNu (STVelOCh/r^lO -10 /^) in different 
galaxies as functions of redshift. The continuous line represents an elliptical galaxy (the 
same as figure 3) and the dotted line indicates a spiral (model M2). The observed local 
type la SN rates in ellipticals (the square) and spirals (the triangle) together with the rate 
measured by Pain et al. (1996) at z ~ 0.4 (the pentagon) are shown. 



